Resonant transmission of light through finite chains of subwavelength holes 
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In this paper we show that the extraordinary optical transmission phenomenon found before in 2D 
hole arrays is already present in a linear chain of subwavelength holes, which can be considered as the 
basic geometrical unit showing this property. In order to study this problem we have developed a new 
theoretical framework, able to analyze the optical properties of finite collections of subwavelength 
apertures and/or dimples (of any shape and placed in arbitrary positions) drilled in a metallic film. 
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After the discovery of extraordinary optical transmis- 
sion (EOT) through 2D square arrays of subwavelength 
holes in an optically thick metallic film Q , several works 
have appeared in order to understand the basics of this 
remarkable phenomenon. From the theoretical side, the 
studies can be divided in those considering the simpler 
ID analog of arrays of subwavelength slits El and 
the 2D arrays of holes |a, la □> lal ■ Several of these works 
explained EOT in terms of the existence of surface elec- 
tromagnetic (EM) resonances, something pointed out by 
the original experiments Q and definitely corroborated 
by recent experiments However, a question that still 
remains open is what is the minimal system that shows 
EOT, which is interesting both from the basic point of 
view and for possible future applications. 

In this letter we move a step forward in this direction 
and consider the optical transmission properties of finite 
chains of subwavelength holes, a basic structure with less 
symmetry than the original 2D array which, up to our 
knowledge, had not been considered before. Here we 
show that EOT is also present in these ID finite systems. 
As an important byproduct, we develop a formalism ca- 
pable of treating the optical properties of even thousands 
of indentations (with any shape and placed arbitrarily) 
in metal films, something not possible with the present 
numerical methods, which are restricted to just a few of 
such indentations. 
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FIG. 1: A chain of circular holes in a metal film with thickness 
h, with a schematic representation of the terms appearing in 
the theoretical formalism presented in the text. 



Let us first present the formalism, which is a non-trivial 
extension of the simpler one developed for sets of ID in- 
dentations and that was successfully applied 0, ^| for 
the understanding of enhanced transmission and beam- 
ing of li ght in single apertures flanked by periodic corru- 
gations |12t Il3| . Here, we analyze the EM transmission 
through a planar metal film (with finite thickness h and 
infinite in the x-y plane) with a set of indentations at 
both input and output interfaces. These indentations 
may be either holes or dimples. Furthermore, each one 
of them may have any desired shape and may be placed 
in any position we wish. The only approximation in the 
formalism is that the metal is treated as a perfect con- 
ductor (e = — oo). Then, within this approximation, the 
method is virtually exact. We have demonstrated in pre- 
vious works that this model captures the basic ingredi- 
ents of the enhanced transmission phenomena, being even 
of semi-quantitative value in the optical regime for good 
conductors like silver or gold [l4|. Additionally, results 
obtained within the perfect conductor approximation are 
scalable to different frequency regimes. 

In our method, we assume a rectangular supercell, with 
lattice parameters L x and L y , along the x and y axes, 
respectively. This supercell may be real (if we are con- 
sidering a bona-fide periodic system) or artificial, if the 
system is finite. In this latter case, the limit L x , L y — > oo 
must be taken. 

For an incident plane wave with parallel wavevector 
fcn|15| and polarization do, the EM field at z = 0~ (at 
the metal interface in which radiation is impinging on) 
can be written, in terms of the reflection coefficients rt , 

' KG ' 

as 

\E{z = Q-) > = \k a >+J2 r kJ^ a> 

ka 

-u z x \H(z = Q-)> = Y^Jkoao > -J^r^Jka > 

ka 

where we have used the Dirac notation, and expressed 
the bi-vectors E = (E x ,E y ) T and H = (H x ,H y ) T 
(T standing for transposition) in terms of the EM 
vacuum eigenmodes, \ka >. The expressions for 
these EM vacuum eigenmodes in real space are: < 
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r\kp >= (k x , k y )' r exp{ikr) / W L x L y \k\ 2 and < r\ks >— 
(—k y , k x ) T exp(ifcr) I ■J L x L y \k\ 2 . The electric and mag- 
netic fields in Eq. (1) are related through the admittances 
Yg s = kz/k^ and = k^/k z (for s- and p- polarization, 
respectively), where fc w = uj/c (u) is the frequency and c 
the speed of light) and \k\ 2 + k 2 = k 2 . Notice that, 
according to Bloch's theorem, k = ko + G, G being a 
(supercell) reciprocal lattice vector. 

In the region of transmission, the electric field at z = 
h + can be expressed as a function of the transmission 
amplitudes t? as \E(z = h + ) >= Y^kaH ""*> f rom 
where the magnetic field can be readily calculated. 

The EM fields inside the indentations can be written, 
in terms of the expansion coefficients A a , B ai as: 

\E(z) > = Y,\ a > [A a e iq "- Z + B a e- lq *" z ] (2) 

a 

-U Z x \H(z) > = ]T |a > Y Q [A a e iq '" z - B a e-^" z ] 

In the previous equations, a runs over all "objects", 
which wc define as any EM eigenmode considered in the 
expansion. An object is, therefore, characterized by the 
indentation it belongs to, by its polarization and by the 
indexes related to the mode spatial dependence. All that 
is required to be known are the electric field bi-vectors 
| a >|l6j and the propagation constants q za associated 
to the objects, as the admittance Y a — q za jk u for TM 
modes, while for TE modes Y a = k u /q za . For inden- 
tations with such simple cross sections as rectangular 
or circular, the required expressions for \a > and q za 
can be found analytically [l7| : otherwise, they can be 
numerically computed [18j. By matching the EM fields 
appropriately on all interfaces, we end up with a set of 
linear equations for the expansion coefficients. We find 
it convenient to define the quantities E a = A a + B a 
and E' a = A a e lq ^ h + B a e~^"' h , which are the modal 
amplitudes of the electric field at the input and output 
interfaces of the indentations, respectively. Objects cor- 
responding to holes are represented by two modal ampli- 
tudes, while objects in dimples require only one (as, in 
this case, A a and B a are related through the boundary 
condition at the dimple closed end). The set {E a ,E' a } 
must satisfy: 

(e a + G aa )E a + 22 GapEp + G^E' a = I a 

(e 7 + G 77 )i? 7 + G 1V E' U + GY f E 1 = (3) 

The different terms in these "tight-binding" like equa- 
tions have a simple interpretation: I a = 2 < k^o^a > 
takes into account the direct initial illumination over ob- 
ject a. e a is related to the bouncing back and forth of the 
EM-fields inside object a and is e a = Y a (l+$ a ) / (l-$ a ) 
where, for holes, (fr a = exp(2iq za h) and the same ex- 
pression applies for dimples, but replacing h by W a , the 



depth of the dimple. The main difference between holes 
and dimples is the presence of G^, which reflects the 
coupling between the two sides of the indentation. For a 
hole G£ = 2Y a $ a /(l - $ a ), while for a dimple G^ = 0. 

The term G a p = J2j: a Yj: a < a\ka >< ka\f3 > con- 
trols the EM-coupling between indentations. It takes 
into account that each point in the object (3 emits EM 
radiation, which is "collected" by the object a. If the 
system is periodic, G a p can be calculated through the 
previous discrete sum, by including enough diffraction 
waves. If the considered supercell is fictitious, the limit 
L x , L y — > oo transforms the previous sum into an integral 
over diffraction modes. It is then convenient to calculate 
G a /3 through G a p =< a|G|/3 >. The integral defining 
the dyadic G(f\\,fu') =< r|G|?Y > can be evaluated, 
obtaining 

Gy (f|, , f{) = g(d) 6a + (2*« - 1) (4) 

where i,j can be either x or y, Sij is Kronecker's delta, 
d = /c w |r|| — f*i/| and g(d) = —ik u exp(td) / (2nd) is propor- 
tional to the scalar free-space Green function associated 
to the Helmholtz equation in 3D. 

G turns out to be the in-plane components of the 
EM Green function dyadic associated to an homogeneous 
medium in three dimensions 01- As a technical note, 
the calculation of G a p for objects in the same indenta- 
tion from < a|G|/3 > suffers from the problems associated 
to the divergence of Gat \f {l - r|h0 [H, E3, EH , and 
we evaluate it directly from its integral over diffraction 
modes. 

Therefore, our method reduces the calculation of EM 
fields into finding the EM field distribution right at the 
indentation openings, which is extremely efficient when 
the openings cover a small fraction of the metal sur- 
face. By projecting these fields into indentation eigen- 
modes, convergence (as a function of number of eigen- 
modes needed) is reached very quickly, especially in the 
subwavelength regime. Notice that, although derived in 
order to treat finite systems, Eq.© can also be used for 
an infinite periodic array of indentations, by imposing 
Bloch's theorem on the set {E a , E' a }, something we will 
make use of when analyzing the case of an infinite linear 
chain (see below). 

Once the self-consistent {E a ,E' a } are found, it is 
straightforward to find all expansion coefficients, and 
from them the EM fields in all space. For the to- 
tal transmittance through subwavelength holes, we find 

We have tested our formalism against published results 
for two extreme systems: a single circular hole [2^] and a 
square array of square holes |7| . In both cases we recover 
the known results, showing that our method is free from 
numerical instabilities and that, in finite systems, there 
are no spurious effects related to the L x , L y — > 00 limit. 
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FIG. 2: (a) Normalized-to-area transmittance (see text) ver- 
sus A/d for a linear array of TV holes with a/d = 0.25 and metal 
thickness h — a. Top inbox: transmittance peak value, as a 
function of N, dashed line showing the value obtained for an 
infinite chain. Bottom inbox: Normalized-to-area transmit- 
tance versus A/d calculated for an infinite 2D hole array with 
the geometrical parameters defining the ID arrays, (b) Gs + e 
and Gv (see text for the definition of these magnitudes) as a 
function of A/d. 



In the rest of this letter we apply this formalism to 
the study of linear chains of subwavelength circular holes 
(see Fig.l). For proof of principle purposes, we choose 
holes with radius a/d — 0.25 perforated in a metallic 
film of thickness h/a — 1, d being the first-neighbor dis- 
tance between holes. These are typical geometrical val- 
ues in experiments in 2D hole arrays. Fig. 2a shows the 
evolution of the transmittance versus wavelength as a 
function of the number of holes, N |2^|. The incident 
plane wave impinges normally, and is polarized with the 
E-field pointing along the direction of the chain. For the 
other polarization, the boost in transmittance is negli- 
gible. In Fig. 2a, the total transmission is normalized to 
the hole area and then divided by the corresponding N . 
The most interesting feature of these spectra is that, as 
N is increased, a transmittance peak emerges at A close 
to d, showing that enhanced transmission is also present 
in linear chains of subwavelength holes. The transmit- 
tance peak value, T max , grows almost linearly with N 
(for small N) eventually reaching saturation (see top in- 
box of Fig. 2a). In order to gain physical insight into 



the origin of this transmission resonance, it is interest- 
ing to analyze a simplified model for the infinitely long 
chain. In this model, only one mode per hole is consid- 
ered: the least evanescent mode with the electric field 
pointing mainly along the chain axis. Therefore, for all 
wavelengths, each hole behaves as a small dipole. In this 
geometry, Bloch's theorem implies E a = Eexp(ik x ad), 
and E' a = E' exp(ik x ad), which renders the system of 
equations (3) easily solvable. For the case we are consid- 
ering (normal incidence, k x = 0) this procedure yields, 

[{e + Gsf -G 2 v ]E = I(e + Gs) (5) 

where G s = G aa + Y*p^ a G af3, G v = G v a , and e = e a . 
Fig. 2b shows that the spectral location of the transmit- 
tance peaks for the infinite chain coincide with the cuts 
between |e + Gs\ and Gy|, implying that the origin of 
EOT relies on a resonant denominator. Therefore, EOT 
is associated, as in the case of 2D hole arrays, to the ex- 
citation of coupled surface EM resonances, which radiate 
into vacuum as they propagate along the surface. This is 
yet another instance of surface EM modes (in this case 
a leaky mode) appearing in a perfect conductor due to 
the presence of an array of indentations |2^. The results 
for finite chains (top inbox of Fig. 2a) show that this EM 
resonance is characterized by a typical length, Lp (for 
the case considered Ld ps 80c?) . When the size of the 
finite chain is smaller than L jj , the resonance is not fully 
developed and the transmittance is smaller than the one 
obtained for an infinite chain. However, for large enough 
N , the system can effectively be considered as infinite and 
its associated transmittance peak approaches the asymp- 
totic value (with a 1/N contribution coming from holes 
located at a distance « Ljj/2 from the chain ends). We 
have found that Ljj is completely governed by the geo- 
metrical parameter a/d, increasing as a/d decreases. 

It is also interesting to compare the peak value calcu- 
lated for an infinite linear array (T^ x w 4, see Fig. 2a) 
with the one obtained for an infinite 2D hole array with 
the same geometrical parameters (T^ x w 5, see bottom 
inbox of Fig. 2a). This is, when going from a linear chain 
to a 2D hole array, the transmittance per hole is only 
increased by 25%. This means that the basic unit of the 
extraordinary transmission phenomenon observed in 2D 
hole arrays is the linear chain of holes, and that the 2D 
hole array can be considered as an array of weakly cou- 
pled ID arrays. Although in this paper we only show 
results for normal incident radiation, we have checked 
that this conclusion remains valid provided the in-plane 
component of the incident E-field points along the direc- 
tion of the chain. 

In order to study how the EOT in 2D hole arrays 
develops from the EOT in linear chains, we have cal- 
culated the transmittance of a collection of finite linear 
chains separated by a distance d. Figure 3 renders the 
normalized-to-area transmittance versus wavelength, for 
different stripes formed by several (ranging from 1 to oo) 
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FIG. 3: Normalized-to-area transmittance versus X/d for 
several stripes formed by M linear chains of 41 holes with 
a/d = 0.25 and h/a = 1. M varies from 1 to 5. Also the 
limiting case of M — oo is presented. Inset: Transmittance 
peak value for the case of two linear chains (each one of them 
with the same parameters as before), as a function of their 
distance, d c . 



chains of 41 holes. As shown in the figure, the stronger 
effect appears when going from 41 x 1 to 41 x 3, suggesting 
that the EM coupling between chains of subwavelength 
holes is very short-ranged. The short-range nature of the 
inter-chain interaction is more clearly seen in the inset to 
Fig. 3, which renders the transmission peak value, T max , 
through two finite linear chains as a function of the dis- 
tance between them, d c . As this inset shows, for this set 
of geometrical parameters, the chains are already prac- 
tically uncoupled when d c = 3d, the maximum coupling 
being for d c rj d. 

When linear chains are added up to the structure, the 
transmittance peak shifts to longer wavelengths and its 
maximum value increases. The increase in transmittance 
is due to the fact that each chain takes advantage from 
the re-illumination coming from other chains. The peak 
redshift is related to the corresponding decrease of fre- 
quency of the stripe surface EM mode, due to a reduc- 
tion of its lateral confinement. More precisely, in this 
case the stripe surface mode that couples to the inci- 
dent plane wave with k = is essentially the sum, with 
equal phases, of the single chain leaky modes. In in- 
finite 2D hole arrays, two EOT peaks originate from 
the resonant coupling (via the holes) of the surface EM 
modes, with a narrow transmission peak appearing close 
to A « d (see bottom inbox of Fig. 2). Notice that, in 
stripes of chains of finite length only one peak is clearly 
resolved, the second expected peak appearing as a shoul- 
der in the transmission curve (see Fig. 3). This suggests 
that finite size effects may prevent the development of 
the narrowest peaks. Interestingly, the addition of linear 
chains also provokes the birth of a minimum in the trans- 
mittance spectrum, appearing at a wavelength slightly 
smaller than d. Eventually, in 2D hole arrays, this mini- 



mum leads to a "Wood's anomaly" , appearing just when 
a propagating Bragg diffracted wave becomes evanescent. 
In this case the reciprocal lattice vectors involved are 
(±l,0)27r/d, and Wood's anomaly appears due to a di- 
vergence in the EM density of states corresponding to 
those wavevectors. Notice that, in a linear chain, as the 
diffracted field contains a continuum of k y components, 
the divergence in the density of states is smeared out; cor- 
respondingly no Wood's anomalies appear in this case. 

To summarize, we have found that extraordinary opti- 
cal transmission phenomenon is already present in a sin- 
gle finite chain of subwavelength holes in metallic films. 
For a chain, the transmittance per hole is comparable 
to that found in 2D arrays; therefore, the single chain 
can be considered as the basic entity of EOT and then 
2D arrays can be seen as a collection of weakly EM cou- 
pled chains. As a byproduct we have developed a new 
theoretical framework which is able to treat the optical 
properties of even thousands of indentations (holes or 
dimples) placed arbitrarily in a metallic film. This nu- 
merical tool will surely help in the design of new arrange- 
ments of subwavelength objects pursuing specific optical 
functionalities. 
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